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THE RADIATION FIELD IS A FOURIER INTEGRAL 

OPERATOR 
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1. Introduction 



In this note, we exhibit exphcitly the form of the "radiation field" of 
F. G. Friedlander on two different types of manifolds: scattering manifolds, 
Ph i and asymptotically hyperbolic manifolds. The former class consists of man- 

■ ifolds with ends that look asymptotically like the large ends of cones, and 
(-H , includes a large class of asymptotically Euclidean spaces, while the latter 

consists of spaces that resemble the hyperbolic space at infinity, and in- 
cludes quotients of hyperbolic space by certain groups of motion. In both 
cases we assume that there are no trapped geodesies. The radiation field is a 
measurement of the (rescaled) asymptotic behavior of solutions to the wave 
equation, viewed from the point of view of a rescaled time coordinate which 
I in the asymptotically Euclidean setting is simply s = t — r, and restricted 

■ to the sphere at infinity. In particular, in M", we define 
O ■ T^+{s, 0, z') = lim r("-^)/2y(^ ^ ^/^ 



' where V{t,z,z') = (cost\/A, sint-v/A/VA) is the solution operator to the 

. wave equation. Friedlander showed in [j^, that 7^+ is in fact a a translation 

I representation of the wave group in the sense introduced by Lax and Phillips 

' ^21- (Sa Barreto I^HOj subsequently showed the unitarity of this map.) 

In this paper, we show that the radiation field 7^+ on a manifold X which 
^ I is either a scattering manifold or an asymptotically hyperbolic manifold, 

has as its Schwartz kernel a Lagrangian distribution, associated to the conic 
Lagrangian defined by the graph of a "sojourn relation" relating points 

o 

in T* X to points in r*(]R x dX), where dX is the boundary at infinity. 

o 

In the simple example where X is just Euclidean space, the graph of the 
sojourn relation maps T*M" to T*(M x S'^~'^) roughly as follows: given 
(z,C) G S'*(M"), let 7 be the unique unit speed geodesic passing through it. 
We map {z,C,) to the base point (s,0) G R x where s is given by the 

"sojourn time" or the limit of t—r along the geodesic, and is the asymptotic 
direction in 5"^^. The fiber variables then measure the angle of contact the 
geodesic makes with S"""^ in a way made precise below. This sojourn relation 
on scattering manifolds was previously investigated by Hassell-Wunsch jH] 
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in the context of the fundamental solution to the Schrodinger equation on 
scattering manifolds; it is also closely related to the sojourn time defined 
by Guillemin in in the study of the high frequency asymptotics of the 
scattering matrix. Note that the sojourn time s is more or less just the 
"Busemann function" used in differential geometry. 

In the special case in which we locally have a finitely many geodesies 7n(i) 



beginning at a point z £X with asymptotic direction 9, the construction is 
simpler: TZ^ is a conormal distribution with respect to the hypersurfaces 
s = Sn{z,0) = limt — r(7„(t)). In this case, Fourier transforming yields a 
simple result about the high-frequency asymptotics of the scattering Poisson 
operator (better known, in the asymptotically hyperbolic the Eisen- 

stein function) . This result is a weak generalization of a result of Guillemin 
0, who proved on compactly supported perturbations of Euclidean space 
that not only do the scattering operator and Poisson operator have the form 
discussed here, but a composition of FIO's gives the asymptotics of the scat- 
tering matrix as well. We are unable to perform this composition owing to 
the local nature of our results in z, the variable in M" (the location of the 
initial pole of the fundamental solution). Our results are also weaker in the 
sense that we obtain only distributional asymptotics of the Poisson opera- 
tor, i.e. we must mollify by convolution with the inverse Fourier transform 
of a compactly supported cutoff in order to describe the asymptotics. On 
the other hand the results presented here are novel insofar as we do permit 
global perturbations of the metric and folded sojourn relations, and we treat 
the asymptotically hyperbolic CQjSG clS well. 

In the case of obstacle, rather than metric, scattering, results relating the 
scattering matrix and the sojourn time were initially obtained by Majda 
[13j . and used in the solution of inverse problems. For further applications 
of sojourn-time methods in inverse obstacle scattering, see the survey by 
Petkov and Stoyanov [TS] . Similar results to ours in the case of semiclassical 
scattering on M" have also been obtained by Robert and Tamura jl9j : more 
recently, Alexandrova P has extended these results to show that even if the 
Hamilton flow is degenerate, the scattering matrix is a semiclassical FIO. 

We now discuss simple examples in which the Poisson operator is explic- 
itly known and the appearance of the sojourn time in the high-frequency 
behavior is clear. 

Example 1. On M"', the kernel of the scattering-theoretic Poisson operator, 
evaluated at (z, 9) gR x 5"~\ is just 



(see ^ni); this is the operator mapping high-frequency incoming scattering 
data to a generalized eigenfunction with eigenvalue (see ^2. 21 for a precise 
definition.). 
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On the other hand, there is a unique geodesic beginning at z with as- 
ymptotic direction : it is just z ^id. The sojourn time along the geodesic 
is 

hm|t - z ^tQ\ = -Q ■ z, 
hence exactly the phase of the adjoint of the Poisson operator. 

Example 2. The kernel of the Poisson operator on hyperbolic space is better 
known as the "Eisenstein function;" it is given by 

Ei^ + i\, y, z) = hm (^^ + i\ x, y, z) 

where R{X) is the resolvent, normalized to be (A — — ^)~^ and where 
we work in the usual coordinates on the half-space with defining function x. 
On H^, the resolvent is just (see equation (6.8) of ^14J: 



R[- + iX) = C- 



■ 2 sinh 5 

where 5 is the hyperbolic distance. The phase is thus —5, which is asymptotic 
to 

x2 + {x'f + {y- y'f 

- log -, , 

xx' 

hence, switching primed and unprimed variables, the phase of £'(^-|-iA, y', 
is just 

(p{x, y, 2/ ) = - log . 

X 

Now we compare this phase to the sojourn time. Given a point z £ 
and y' G OH^, there is a unique geodesic 7(t) starting at z, and approaching 
y' as i — > oo (this is the analogue of the "asymptotic direction" 6 G S^~^ in 
the Euclidean case). We now define the sojourn time for this geodesic as 

S = S{z, y') = limt -|- logx(7(t)) 

Using translation invariance in the boundary variables, it suffices to compute 
with y' = 0. A unit speed geodesic from an arbitrary point {x, y) S to 
(0, 0) G dM? can be parametrized as 

X = Csech(f + to), y = Ctanh(t + to) + D, 

whence we compute (p{x, y, 0) = — log 2 — log C + to. On the other hand, the 
sojourn time along such a geodesic is 

S{x,y,0)= lim t log x(t) = log 2 log C - to, 
hence agrees with minus the phase (j){x,y,0). 
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2. Scattering manifolds 

2.1. Radiation field. Let X he a C°° compact manifold with a boundary. 
A scattering metric (defined originally by Melrose is a metric on a 

manifold with boundary that can be brought to the form 

^ , dx'^ h 

with X a boundary defining function and h a smooth tensor that restricts 
to X = to give a metric ho on dX. This form is modeled after the metric 
on asymptotically Euclidean space, radially compactified. It was shown by 
Joshi-Sa Barreto [HI that there is a normal form for a scattering metric: 
in which a neighborhood of the boundary admits a product decomposition 
[0, e) X dX with local coordinates {x, y) in which g takes the form (|2.1|) 
with h = h{x,y,dy) a smooth family in x of metrics on dX. Note that in 
these product coordinates, the rays y = constant are (infinitely extended) 
geodesies. We henceforth assume that our metric is in normal form. We 
further make the geometric assumption that there are no trapped rays in 

o 

X . 

Following Friedlander we define the forward radiation field on a scat- 
tering manifolds by 

7^+(/l,/2)(s,y) = {x-^''-'^/^DtHu){s + x-\x,y)\,^o 

where H = H{t) is the Heaviside function and u = u{t, x, y) is the solution 
to 

(2.2) Uu = ^, (tx,An)rt=o = (/i,/2). 

Friedlander showed that 7^+(/i,/2) G C°°(M x dX) provided /i,/2 are 

o 

smooth and compactly supported in X ■ Thus the Schwartz kernel 

7^+(s,y,z) 

o 

is defined on M x dX; here we have used coordinates (s, y, z) on M x dXx X ■ 
First we will show that the radiation field is an FIO and then, under some 
nondegeneracy conditions, we will compute its symbol. We use coordinates 

o 

on T*(M X dXy. X) defined by the canonical one-form 

sda + ydrj + z dC, 
and will employ the notation C = C/ICI- 

Theorem 2.1. Let X be a nontrapping scattering manifold. Then 
1Z+ G (/^(R X dXx X, A'; 1^1/2), /°(R x dXx X, A'; l^^/^)), 
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where fl is the density bundle onMxdXx X, and A'^ is the conic Lagrangian 
manifold associated to the graph of the "sojourn relation": 

A± = A+ U 

A± = i 2, C, s = lim t- x"^(exp^(tC)), cr = ±\C\, 

y= lim 2/(exp^(tC)),r/j= lim ±\C\{ho)ij{dy^ /dx){exp-.{tC))\. 

Proof. We know from the non-trapping assumption and from Theorem 1.1 
of P] that the kernel of cos t\/A, which is the solution to (|2.2j) with /i = 5{z) 
and /2 = 0, a sum of forward and backward parts (corresponding to a choice 
of lb in our notation below), each of which is a Lagrangian distribution in 

/"i (M X X X X; C±) , where 

C± = {(t, z, z', r, C, C) : C), (^', C) e \ 0, (t, r) G T*M \ 0, 

T = ±^/L{^), {z, C) = ^^\z', C')}. 

Here L{z,(^) denotes the symbol of A and <I>*(z',C') denotes the flow along 
Hl in T*X \ 0. We want to understand what happens to the kernel of 
cos ty/A if one first makes the change s = t — 1/x and then takes the limit 
X ^ 0. So we need to understand the effect of these operations on C± and 
the distribution associated with it. 

As in let P = ^"^"'^"'"^^/^□x^""^)/^ and change variables, replacing 
t by s = i — 1/x. We find as in that 

1 /2 

with ^(x, y) = dx log \ h\ ' and A/j is the nonnegative Laplacian on dX with 
respect to the metric h. The symbol of P in these coordinates is given by 

p = -2^a - x^i^ - h{x, y, ry) 

and the Hamilton vector field by 

((zjC), the coordinates in the right factor, are left invariant by the flow). 

If Vo = X ~Uo, the equation 1)2. 2|) with initial data fi = 6{z) and /2 = 
becomes 

n — l 

(2.4) PVo = 0, Vo = x- — 6{z), DtVo =0, x > 

X x 

We remark that we are ignoring half-density factors, as they are irrelevant to 
this construction. Notice that the operator P extends to x < as a strictly 
hyperbolic differential operator P and that s = — ^ is a space-like surface for 
P. One can think of this as being an extension of P to the double manifold 
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X"^ = {XUX)/dX. Therefore it follows from the existence of a fundamental 
solution to the Cauchy problem for strictly hyperbolic operators, see for 
example Theorem 5.1.6 of [2], that the solution to (|2.4|) . with P replaced by 
P, is a Lagrangian distribution Vq of class 

(2.5) Vo e (rx Xx X;C 

where C denotes the Lagrangian in T*(M x Xx X) obtained by flowing 

N*{z = z', t = 0} n Sp C r*(IR X Xx X) along the integral curves of 
Hp, where p is the principal symbol of P, and Sp denotes the characteristic 
variety of P. By the uniqueness of solutions to the Cauchy problem, the 
restriction of Vq solution to x > is equal to Vq, the solution to (|2.4jl . 

We remark that the extension P, and consequently the definition of £ in 
{x < 0}, are not unique. However the extensions to {x = 0} are. 

Since N*{{x = 0}) H C = 0, U = VqI^^q is a Lagrangian distribution of 

o 

order 1 on M x dXx X with respect to CI^^q. Hence for any compactly 
supported distribution u, 7^+(m, 0) = DtUu where C/ is a Lagrangian of 
order with respect to C\^^q. 

To identify C\^^q geometrically we observe that under the flowout of Hp, 
a is conserved, hence by homogeneity we need only consider a = ±1 (note 
that o" / on £ \ 0). Let li± denote the two components corresponding 
to different signs of a. By definition, we certainly have s = \\m.i^±ao^ — 
x^^), and y is the limiting location of geodesic flow (forward or backward 
according to sgncj) in dX. Furthermore, setting a = ±1 we have dy'^ /dx = 
{-20r]j)/{-2{±l + x'^i)), which approaches ±h^^ rjj as X ^ 0. Thus the 
restriction of C± equals A± as defined in the statement of the theorem. 

Similarly, solving 

(2.6) PUi = 0, C/irt=o = 0, DtUo\t=o = S{z) 

gives a Lagrangian of order —1 when restricted to the boundary. Hence for 
any distributions /o , /i , 

^+(/o,/2) = (At/o/o, A?7i/i)r,.=o 
is a FIO of the asserted kind. □ 

It remains to calculate the symbol of 7^+ and for this we need some extra 
assumptions. So we suppose further that for all z contained in an open set 

o 

Ui in X, all y contained in an open set U2 in dX, there exist a finite number 
of unit speed geodesies 7,1(2:, y, i), n = 1, . . . , such that 

7n{z,y,0) = z, lim 7„(z,y,t) = y, 

1^+00 

and such that the transformation 
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mapping 

T{dX) ^ T{S{X)) 
is invertible for all n, and y £ U2, z £ Ui. Subject to these assumptions, the 
Lagrangian A is projectable onto the z, y variables. Letting 

Sn{z,y) = Jim. t - x~^{jn{z,y,t)) 

I— > + oo 

(the "sojourn times"), we then find that 7^+(0, 6z) is conormal to the surfaces 
s = Sn{z,y) in Rs x Ui x U2. 

Theorem 2.2. Subject to the nondegeneracy assumptions above, the symbol 
ofTZ-^-, evaluated at (z,y,a) G N*{s = 5„} equals 



(2.7) -i^ 



^1/2 



I ds dhy dg^ 



1 1/2 . / (n+l)/2 Jn-l)/2^ 



dc 

where kn is the number of conjugate points encountered by ^n{z,y,t) with 



t G (0,00) and dy^''\z,y)/dC 



is short for the Jacobian of the map 



SliX) ^ dX 

given by the limit of geodesic flow, evaluated at the initial codirection ofjn- 

Note that we have written the Jacobian factor in the above form to empha- 
size the analogy with the differential scattering cross-section in Guillemin's 
results [Z]. 

o 

Proof. Let OG C S*{X) denote the outgoing set, i.e. the set on which 
dx/dt < along the bicharacteristic flow. Owing to our nontrapping as- 

o 

sumption, the cosphere bundle of a compact set S*{K) C S*{X) eventually 
maps into OG under the bicharacteristic flow after time T ^> 0. Let us fix 
such K and T, with K chosen such that our nondegeneracy assumption 
holds for all geodesies beginning in K with limits some open set in y. 

Net Z be a compactly supported pseudodifferential operator with WF'Z C 
OG, WF'{1 - Z)n OG' = 0, with OG' a conic subset of OG with compact 
projection, chosen so that the flowout of S*K for time T lies inside OG'. 

Let Zj be a microlocal partition of unity over K such that for all j, and all 

z £X, WF' Zjrnr~^{z) contains at most a single point such that the geodesic 
emanating from this point ends up at y G dX. By our hypotheses on the 
nondegeneracy of geodesies, we may further choose Z supported sufficiently 
close to y that if {z' , (^') and {z, (") are canonical coordinates on the cotangent 
bundles of left and right factors, then {(^', z) are coordinates on 

A± n TT^WF'Z n vrljWF'Zj. 

We may further arrange, by working sufficiently close to dX, that there are 
no conjugate points for bicharacteristics beginning on WF'Z : examination 
of the Hamilton flow in "scattering coordinates" shows that the tangent 
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vectors to geodesies emanating from WF'Z approach —x'^dx, and while the 
sectional curvature of a scattering manifold is O(x^), the sectional curvature 
of a plane containing x^dx is O(x^). Hence certainly the sectional curvature 
along a plane containing the tangent to the geodesic is o(x^). By a simple 
variant on the Rauch comparison theorem, this is sufficient to ensure the 
absence of conjugate points. (See ^Sl for a description of the bicharacteristic 
flow in scattering coordinates and, for instance Theorem 4.5.1 of for an 
account of the relevant comparison theorem.) 

Now let W{t) denote the propagator for the wave group. The symbol of 
W{T)Z'- is 



ici; 



1 



where kn is the number of conjugate points encountered. Changing to coor- 
dinates (C', z) gives 

'1/2 



1 



dc 



dC,' dz 



1/2 



1 

IC'I, 



IC1 



1 



By invariance under the bicharacteristic flow (see for instance Proposition 
4.3.1 of (2), the symbol of U+Z is 



1 



d{s,y) 



-1/2 



ds dy dz'\^^'^ {\a\, 1). 



We split the coordinate into |C|')C' (where the latter should really be 
regarded as n — 1 components of C,'). By homogeneity of the flow we have 
9s/5|C| = 1, while dy/d\C,\' = 0, hence the symbol of IZj^-Z can in fact be 
written 



d{y) 



-1/2 



1 /2 

\dsdydz'\ {\a\,l). 



The difference 

■R+{s)Z'j - n+{s - T)ZW{T)Z'j 

is a smoothing operator. Hence applying the calculus of FIO's to the above 
results and patching together the partition of unity Zj, we find that the 
symbol of 7^+ is given globally by 1)2. 7() . □ 

2.2. Poisson operator. Let P{X) denote (the Schwartz kernel of) the Pois- 
son operator, i.e. the operator such that for any g G C°°{dX) there exists 
u G C°°{X) with (A - X^)u = and 

with g' G C°^{dX), u' G L'^{X;g). 

As a corollary of the results in the preceding section, we conclude the 
following: 
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Proposition 2.3. Suppose as in the preceding proposition that there exist 
a finite number of nondegenerate geodesies 7n from z to y. Let (p G 
Then as |A| oo, 

^ / A \ - -1/2 

<i){X) * P(A)*(y, z) ~ <^(A) * Yl i''"e'^^- i A I 



n=l 



27ri y 



dC 



Proof. We simply use the fact that 

P(Ar(y,z) = -2^7^+(0,<5, 

(see in])- 



□ 



Note: if we knew more about energy decay, along the lines of having good 
estimates for the energy norm 

d voL 



dV{t)f 



dt 



E{tA) = \\V{t)f\\H^^,_^^= [ l\VVit)f\' + 

with initial data / compactly supported, we would be able to get better 
estimates for the decay as s ^ +oo of TZ+, and hence drop the mollifier (f> 
from the statement of this proposition. Friedlander proves in 4 that 

lim lim E{t, A) = 

A— »oo t— >oo 

where / is finite energy initial data; we would need a good deal more how- 
ever. 



3. Asymptotically hyperbolic manifolds 

An asymptotically hyperbolic, or conformally compact manifold is a mani- 
fold X with boundary equipped with a metric g and defining function x >0 
such that x'^g = H is a smooth metric on X, nondegenerate at dX. We 
further take = 1 on dX, which ensures that the sectional curvatures 

approach —1 at dX (see, for instance, |14j.) 

Notice that H\qx is only defined by g modulo a conformal factor. It is 
shown in see also 0, that given a conformal representative of H\qx 
there exists a unique boundary defining function x such that 

dx^ + h{x,y,dy) . . , 

x^ 

From now on fix these coordinates. The definition of the radiation fields will 
depend on this choice of x, or the conformal representative of -fffg^. 

We refer the reader to for a discussion of the Eisenstein function on 
an asymptotically hyperbolic manifold. 

In this setting, we can prove more or less the same results as in the case 
of scattering manifolds; the analogous results are as follows: 
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Theorem 3.1. Let X be a nontrapping asymptotically hyperbolic manifold 
of dimension n. 
(1) We have 

TZ+ e (/^(M X dXx X,A';l^^/2)^ jO^j^ ^ X, A'; ^^^/^)), 
where Aj_ is the conic Lagrangian manifold associated to the graph of 
(3.1) A± = A+UA_, 



^± = <,z,C,s = lim t + logx{exp^{tC)),a = ±\C\, 

I t^±oo 

y= lim y(exp^(tC)),r/i= lim ±x"^|C|(/io)ii(c?y'/f^2;)(exp^(tC)) 

(^^j Subject to the nondegeneracy assumptions of '>\2.1\. the symbol of TZ^ is 
given by ()2.7() . 

(3) Let + i\,y, z) denote the (transpose of the) Eisenstein function. Let 
(j) G C^(]R). Subject to the nondegeneracy assumptions of '<\2.1\. as |A| — > oo, 



N . , x{n-l)/2 



<^(A) * + ^A, y, ~ 0(A) * ^ .^"e^^^" 

n=l 



5C 

Proof. It is shown in Sa Barreto |2nj if X has dimension n, and u = cos ^tV A 

than the rescaled fundamental solution, v = x~*-"'~^^/^u(s — logx,x,y), sat- 
isfies 



(3.2) 



^ , d { d d \ ^ , d , d n— 1,, 
Pv=[ — [2— + X— \-x/^h + A— + Ax— + ^—A \v = Q 
ox \ OS Ox j OS Ox 1 

dv 



,-(-i)/25(^), _^^^j^^^ = o, x>0. 
The symbol of P is thus — (2^cr + x^^ + xh(x, y, ry)), hence the flow is 



'^U = logX~'^ ■ " 'J'y^lj Q^\s = \ogX 



dh 

Hp = -2{a + xOdcc + if + h + x^)d^ - 2^9, - xHh. 

As in the scattering case the operator P is strictly hyperbolic in x > and 
{s = logx} is space-like. Moreover P has an extension P to a neighborhood 
of {x = 0}. However, P is not strictly hyperbolic at x = 0. So in principle 
we can only guarantee that v satisfies (|2.5j) . where C is defined with respect 
to P in when x > 0. 

Since P can be chosen to be of real principle type, the parametrix con- 
struction of the Cauchy problem can be carried though across {x = 0}. This 
guarantees that in a neighborhood W of dX there exists 

Vq G L--i flR xWx X;C 
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such that 

PVo = f e C°°(M X W) 
(3.3) yo\s=iogx - H^) = 9 G c^d), 

dVo ^ o 

\s=iogx = h € C^{X)- 

Moreover by finite speed of propagation one has that / is supported in 
s > sq, for some sq. Again using the hyperbohcity of P in x > 0, there exists 
V eC°^{Wn{x> 0}) satisfying 

PV = f in x>0 
dV 

with V supported in s > si. Since / is supported in s > sq and is smooth 
up to dX, and the initial data is compactly supported, the proof of Theo- 
rem 2.1 of j2n| (in particular, the extension of the energy estimates to the 
inhomogeneous equation) shows that V has a smooth extension up to dX. 

Therefore the solution to (|3.2() satisfies v <^ [m.'k Xx X; up to dX. 

The transversality of the flow to {x = 0} and the facts that u 7^ 
on the characteristic variety for a; > and Hpa = imply that f fa.=o is 

o 

a Lagrangian distribution of order on M x dXx X with respect to the 
Lagrangian £f^=0) where £ is the flowout of the hft of A^*A C T*{X x X) 
to Sp. We find in this setting that lim^_^o x^^dy^/dx = rjjh^Q /a, hence C \x=q 
has the form ()3.1() . 

The remainder of the proof is the same as in the scattering case, using 
the additional fact from [10] that E{n/2 + i\y, z) = -{i/\)Tn+{Q,5^). □ 
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